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Notes to Table 1. 

(1) RC3, (2) TuUy 1989, (3) P2, (4) UGC, (5) PI, (6) TuUy-Fisher, 
see Kenney et al. 1993, (7) Sage 1993b, (8) Combes et al. 1992, (9) 
G-B, (10) PSC, (11) Sage 1993, (12) Lonsdale et al. 1985. 
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Table 2. OVRO Interferometer Results for NGC 4314 



Parameter 


Value 




Phase center 






RA 
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DEC 
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CO intensity peak 
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DEC 
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CO kinematic center 






RA 
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(J2000) 


DEC 


+29°53'42'.'6 




Inclination 


21°± 20° (adopted) 




p. a. of line of nodes 


115°± 10° 




Field of view (HPBW of primaries) 


65" 




Baselines 


4-77 k A 




Velocity of CO kinematic center (LSR) (Vsys) 


983 ± 5 km s'^ 




Velocity resolution 


13 km s~^ 




Map-specific Parameters 
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Size of synthesized beam 


2'.'29 X 2'.'17 


3'.' 12 X 2'.'7 


Linear size of synthesized beam 
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151 X 131 pc 
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5.63 


9.6 arcsec^ 


Position angle of synthesized beam 


-16.6 


87.7° 


RMS noise in channel maps 


0.028 


0.022 Jy beam-i 


Equivalent TB for 1 Jy beam~^ 


19.4 


11.4 K 


Peak Brightness Temperature (Tb) 


4.9 


3.7 K 


Interferometer flux 


245 Jy km s^^ 




CO flux (r < 8") 


190 Jy km s"^ 




L(CO) (r < 8") 


4.7 X 10"^ K km s'^pc'^ 




Mass of molecular gas (r < 8") 


2.1 X 10*A//q 


2.1 X 10*^ M0 
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ABSTRACT 



NGC 4314 is an early-type barred galaxy containing a nuclear ring of 
recent star formation. We present CO (1-0) interferometer data of the bar and 
circumnuclear region with 2.3" x 2.2" spatial resolution and 13 km s~^ velocity 
resolution acquired at the Owens Valley Radio Observatory . These data reveal 
a clumpy circumnuclear ring of molecular gas. We also find a peak of CO inside 
the ring within 2" of the optical center that is not associated with massive star 
formation. We construct a rotation curve from these CO kinematic data and 
the mass model of Combes et al. (1992). Using this rotation curve, we have 
identified the location of orbital resonances in the galaxy. Assuming that the 
bar ends at corotation, the circumnuclear ring of star formation lies between two 
Inner Lindblad Resonances, while the nuclear stellar bar ends near the IILR. 
Deviations from circular motion are detected just beyond the CO and Ha ring, 
where the dust lanes along the leading edge of the bar intersect the nuclear ring. 
These non-circular motions along the minor axis correspond to radially inward 
streaming motions at speeds of 20 — 90 km s~^ and clearly show inflowing gas 
feeding an ILR ring. There are bright H ll regions near the ends of this inflow 
region, perhaps indicating triggering of star formation by the inflow. 
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1. Introduction 

Inner Lindblad resonances (ILRs) are thought to strongly influence gas motions and 
star formation in the central regions of galaxies. Gas surface densities often peak near the 
ILR (Combes et al. 1992; Kenney et al. 1992), lending support to the idea that the inward 
flow of gas along galaxian bars slows down and piles up between the OILR and the IILR 
(Combes 1988; Shlosman et al. 1989). Studies of the CO (1-0) distribution within the 
central regions of barred spirals (Kenney et al. 1992; Kenney 1995; Regan, Vogel, & Teuben 
1995; Rand 1995; Sakamoto et al. 1995) show a variety of molecular gas morphologies 
near the ILR, including rings, spiral arms, and "twin peaks" with two maxima where the 
dust lanes along the primary bar cross the ILR. In the case of the starburst galaxy NGC 
3504, the CO peaks near the nucleus, probably inside the ILR (Kenney et al. 1993). The 
dynamical and evolutionary differences responsible for this variety are not yet understood, 
but are critical for understanding starbursts and the formation and evolution of compact 
circumnuclear stellar disks (Kormendy 1993) and stellar bulges (Pfenniger & Norman 1990). 

At present, only a handful of barred spirals have been mapped in CO at sufficient 
resolution to study the effects of ILRs on the gas distribution and kinematics. The barred 
spiral galaxy NGC 4314 is a particularly good galaxy to study because it exhibits evidence 
for many features associated with resonances and it is nearly face-on {i ~ 23°). It has a 
large-scale stellar bar of diameter 130" (6.3 kpc) and a prominent circumnuclear ring of star 
formation of diameter 10" (500 pc) that is visible in Ha (Pogge 1989), radio continuum 
(Garcia-Barreto et al. 1991; hereafter G-B), and optical color maps (Benedict et al. 1992; 
hereafter PI). CO (1-0) mapping of NGC 4314 (Combes et al. 1992) at 5" resolution 
revealed the presence of a molecular ring slightly smaller than that of the ring of star 
formation. Outside of these rings, a blue elliptical feature of diameter 20-25" (1 kpc) is seen 
in optical color maps, which may correspond to a ring of relatively young but non-ionizing 
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stars (PI). This feature, which can be seen on the unsharp masked optical image shown in 
Fig. 1 (from PI), is elongated perpendicular to the primary bar, suggesting that it extends 
outward to the IILR. Combes et al. (1992) suggest that in this galaxy star formation is 
propagating inwards. This is in contrast with the evolutionary scenario proposed by Kenney 
et al. (1993) for strong starbursts, in which star formation devours gas most rapidly in the 
center, ultimately forming a ring of gas near the ILR. 

Inside the H ll region ring, a nuclear stellar bar of diameter 8" (400 pc) is seen in Hubble 
Space Telescope (HST) I band data (Benedict et al. 1993; hereafter P2). This nuclear bar is 
aligned parallel with the large scale bar, which suggests that the nuclear bar has the same 
pattern speed as the primary bar, and lies within the IILR (Binney & Tremaine 1987). 

In this paper, we present new CO (1-0) maps of NGC 4314, obtained at twice the 
spatial resolution as the Combes et al. (1992) data. The observations are described in 
Section 2, while the CO distribution and kinematics are discussed in Sections 3 and 4, 
respectively. In Section 5 we compare NGC 4314 with other galaxies and discuss star 
formation in this galaxy. Some general properties of NGC 4314 are tabulated in Table 1. 



2. Observations 



We observed NGC 4314 in the CO (1-0) line using the Owens Valley Radio Observatory 
(OVRO) millimeter-wave interferometer (Padin et al. 1991) in five configurations between 
April and June 1991. Projected basehnes ranged in length from 10 to 200 meters. At that 
time OVRO consisted of a three element array with 10.4m dishes (HPBW = 65" and SIS 
receivers. The 32 5 MHz receivers give 13 km s~^ resolution and an instantaneous bandpass 
of 416 km s~^. The quasar 1219 -|- 285 was used for phase calibration, Uranus and Neptune 
for flux calibration. 
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After calibration, channel maps were made with both uniform and natural weightings 
using AIPS. Uniform weighting produces higher resolution maps, with less sensitivity, while 
natural weighting provides greater sensitivity with less resolution. The latter reduction 
mode provides higher signal to noise information on the regions in the primary bar of NGC 
4314 containing fainter CO intensity. 

All maps were made as 256 x 256 arrays with 0'/5 pixel"^. Natural weighting produced 
a resolution of 3.1 x 2'.'7 with channel noise 0.022 Jy beam~^. These maps were CLEANed 
to a 0.045 Jy beam~^ level. Using uniform weighting, wc obtained a 2.3 x 2'.'2 beam with 
channel noise 0.028 Jy beam~^, CLEANed to a 0.07 Jy beam~^ level. For both weightings 
we have made primary beam corrections. Fig. 2 presents our uniformly weighted channel 
maps. These maps have absolute positional accuracy of about 0'.'5. Table 2 summarizes the 
interferometer-related parameters . 

The total CO flux observed by OVRO is 245 Jy km s"^ (Table 2), while the CO 
luminosity is 4.7 x 10'' K km s~^ pc^. This luminosity is consistent with the single dish 
observations of G-B and with the Nobeyama Millimeter Array results of Combes et al. 
(1992) (see Table 1). Therefore we have recovered essentially all of the CO flux from this 
galaxy with our OVRO observations. 

3. Morphological Results 

3.1. CO Morphology 

Uniformly weighted and naturally weighted total intensity maps were constructed using 
the AIPS moment routines. Fig. 3 presents the uniformly weighted CO intensity map. The 
strongest CO emission arises from an incomplete, clumpy ring with an off-center minimum. 
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There are five clumps of emission along the ring, separated by gaps. The largest of these 
gaps is at p. a. = 239°. A secondary gap occurs at p. a. = 330°. In addition to the clumps in 
the ring, we find strong CO emission within 2" of the optical center. 

Two symmetrically located features along p. a. ~20° extend out from the ring. These 
are associated with the dust lanes that lie along the leading edge of the primary stellar bar. 
These lanes are detected in CO as they curve around to intersect the nuclear ring. We refer 
to these features as inflow spurs, since the CO velocity field shows clear evidence for inflow 
motions here (§4.2). 

Fig. 4 shows the naturally weighted CO intensity distribution, covering about 4 times 
the area on the sky as Fig. 3. Note the CO clumps extending on either side out from the 
ring along p. a. ~148°. These extensions lie along the leading edge of the primary stellar bar 
(PI). 



3.2. Compctrison of CO to Optical and Radio Continuum 

3.2.1. Maps 

In Fig. 5a we overlay the uniform-weight CO map contours on the optical unsharp 
masked image shown in Fig. 1. The dust and CO are relatively coincident for the CO arcs 
to the SE and NW, while the dust lanes lie outside the CO ring to the NE and SW. Along 
the primary bar, a dusty area is seen to the SE (see Fig. 1), a region called the 'dust bowl' 
in PI. A comparison of optical and CO intensity from the naturally- weighted data for a 
larger area of the galaxy (Fig. 5b) shows weak (~ 3(t) CO sources coincident with the dust 
bowl. Note in both Fig. 5a and 5b the anticoincidence of CO with the outer ellipse of 
relatively young stars. 
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In Fig. 6 the Pogge (1989) Ho; + [N II] map is overlaid on a uniform weight CO 
intensity contour plot. The Wakamatsu & Nishida (1980) Ha-bright knots A (p. a. ~ 90°) 
and B (p. a. ~ 178°) are the second and third strongest Ha sources in the Pogge map. The 
Ha peaks tend to lie outside of the CO ring (see also the profiles in Fig. 7-9). At the 
largest CO gap, Ha is also weak. 

The Ha + [N II] map resembles the 6 cm radio continuum map shown in G-B and our 
observations confirm the results of Combes et al. (1992) that the ring of star formation has 
a larger radius than the molecular ring. 



3.2.2. Profiles 



The profiles in Fig. 7-10 offer a more detailed and quantitative comparison between 
the distributions of CO intensity, blue light {fi^, Pl)j /U{Ha+[N//]) (uncalibrated Ha-|-[N II] 
surface intensity, Pogge 1989), and /ie-H (respectively identified as the best tracer of 
dust and of new stars in PI). Fig. 7, 8, and 9 are primary bar major axis, primary bar 
minor axis, and E-W cuts across the galaxy, respectively. Fig. 10 shows a cut through the 
dust bowl, perpendicular to the primary bar. For these plots the uniformly- weighted CO 
intensity has been converted to (7h2(Mq pc~^) assuming (Bloeman et al. 1986, Kenney et 
al. 1992) 

(Th2(Mqpc^^) = 470 Ico (Jykms~"^arcsec~^) 

Ha is a tracer of the local star formation rate (SFR). These plots also contain the log of 
the ratio of (Ha+[N II]) /CO. In Hll regions, where the observed Ha-|-[N II] emission is 
dominated by Ha and the ionization is due to OB stars, this ratio is proportional to the gas 
depletion timescale, or Star Formation Efficiency. 
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In Fig. 7, the dust lanes ai r — ±5'.'5, seen as reddening in the profile, are very 
prominent in CO. There is also considerable CO emission near the nucleus; Fig. 7, 8, and 9 
show the CO surface density near the nucleus as high as in the ring of star formation. The 
A*(Ha+[N7/]) profiles corroborate that the H ll regions in the ring lie slightly outside the CO 
ring. Dips in the //b-h profile are coincident with the new stars at the Ha ring, r = ±7" . 
The (Ha!+[N II]) /CO ratio shows a relatively higher SFE near the nuclear ring, especially 
to the NW. 

There is also a secondary source in //(Ha+[Nii])) peaking ~ 1" S of the center. This is 
likely due to strong [N II] emission contaminating the Hq;+[N II] image near the center. Keel 
(1983) identified NGC 4314 as having an "[N II]"-type optical spectrum. Measured values 
of [N II]/Ha for the NGC 4314 nucleus range between 1 and 2.5 (Keel 1983; Wakamatsu & 
Nishida 1980; Stauffer 1982; Smith et al. 1987) and Keel (1983) identified NGC 4314 as an 
"[N II] "-type galaxy. The ratios of Ha, R(3, [N II] A6584, [O HI] A5007, and [O I] A6300 
given in Smith et al. (1987) are consistent with the LINER definition given in Heckman 
(1980) & Dahari (1985). Therefore the nuclear peaks in Ha+[N II] and (Ha+[N II])/CO 
may not be due to star formation. 

Note the coincidence of an increase in profile, a decrease in /xb, and a peak in the 
(Th2 profile to the SE. These all argue that CO and dust are associated. How good is this 
correlation? Is there CO only where there is dust? The peak in the profile to the SE 
(Fig. 7) is well-described by a gaussian centered at r = — 5'.'4, with FWHM = 4'/14. From 
the HST data presented in P2, we find the dust lane can also be described by a gaussian 
centered at r = — 5'.'5, with FWHM = 1'.'60. Along this p. a. the OVRO synthesized beam 
can be considered to be a gaussian with FWHM = 2'.'20. Convolving the OVRO beam with 
the dust profile produces a distribution with FWHM = 2'.'70. Evidently, the CO is not 
confined to the dust lane, but spills over in a symmetric distribution. We next deconvolve 
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the OVRO beam from the Fig. 7 (Thj profile and determine that the actual CO distribution 
has FWHM = 3'/51 and a peak = 2550 Mq pc"^ at the dust lane center. 

Along a shce perpendicular to the primary stellar bar (Fig. 8), we again find strong 
CO near the nucleus and near the dust lanes associated with the ring. (Hq;+[N II]) peaks 
outside the CO to the NE, but inside to the SW. The (Hq;+[N II]) /CO ratio picks up near 
r = —4" at a gap in the CO ring. There is an increase in the SFR at r = — 1" , due to the 
observed CO deficiency, or to [N II] contamination from the LINER phenomenon. To the 
SW, the redder dust lane at r = 6" is present in CO as well. The CO maximum at r ~ 4"5 
to the NE is not as pronounced in the color indices as that to the SW, to be expected if the 
SW side of the galaxy is nearer. 

For completeness we provide a set of profiles (Fig. 9) passing East- West through knot 
A, the second most intense Hll region seen in the Pogge (1989) data. The dust to CO 
correlation is not as striking, primarily because dust signatures are weak in and /xb-h- 

There is one other location where comparisons with optical data might prove useful; at 
the "dust bowl" in the primary stellar bar to the SE at r ~ 23" . In Fig. 10 we plot profiles 
passing through the dust bowl, perpendicular to the primary stellar bar. The profiles are 
Hb, A*v-i, and /xb-v (PI) along with a-H.2 derived from the natural weight map (Fig. 3). 
Again, dust and CO are well correlated, with reddening in //y-i matching increases in au^ 
all along the profile. 

We test the reahty of this dust and CO correlation in Fig. 11, showing the relationship 
between the measured CO intensity and Ay across the dust bowl in the range 1 < x < 7" 
. The Ay values were derived in PI from B-V surface color indices. The correlation 
demonstrates the quality of the CO data at the dust bowl location. This suggests that most 
of the weak CO features detected along the primary stellar bar (Fig. 3) are probably real. 
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As a final check of the reahty of this feature in CO, we derive an Np to Ay relationship. 
Prom Fig. 11 we determine 

Ico(Jy arcsec^^ km s~^) = a + b Ay (mag) 

where a = 0.080 ± 0.023 and b = 1.955 ± 0.178. Converting the slope to K km s~^, using 
(Table 2) 1 Jy arcsec"^ = 1.19 K, we find a slope of 2.54 K km s~^ mag~^. Prom Bloemen 
et al. (1986) 

N(H2)(cm-2) = 2.8 x 10^° Ico (Kkms"^) 

Hence, 

N(H2) = 7.1 X 10^° Ay 
and, introducing a factor of two in converting from H2 to Np, 

Np = 1.4 X 10^^ Ay 

in general agreement with a value 

Np = 1.8 X 10^^ Ay 

derived from Savage & Mathis (1979), assuming Ay/E(B-V)= 3.1. The agreement suggests 
that, at least in this part of NGC 4314, a dust screen description is valid (Witt et al. 1992). 

Unfortunately, as discussed in PI, the dust bowl is the only region in NGC 4314 for 
which we can obtain reasonable values for Ay. This is possible only by exploiting the 
symmetry and the relatively uniform stellar population of the primary stellar bar. In the 
more complex circumnuclear region, B-V color variations are due to spatial variations in 
the stellar populations, as well as dust extinction. 



4. Kinematical Results 



4.1. The CO Velocity Field and the Rotation Curve 



In Fig. 12 we show the naturally weighted velocity field, which is consistent with 
predominately circular motion for r < 10" . Prom these data, we have determined the 
systemic velocity Vgys, the dynamical center of the galaxy {a,S), the position angle of 
the line of nodes (p, the inclination i, and a rotation curve V(R). We used the iterative 
technique of Puche, Carignan, and Wainscoat (1991), which employs a tilted ring model. 
Our data was binned into seven rings, each I'.'S in width. Velocity values within each ring 
were weighted by the cosine of the azimuthal angle (relative to the line of nodes) in the 
plane of the galaxy. We first held and i constant and fit for a, 5, and Vgys- Next, we 
held a, 5, and Vgys constant and fit for 0, i, and V(r), fitting the receding and approaching 
halves of the velocity field separately. We included only data within 45° of the major axis. 
Our results are listed in Table 2. 

The inchnation derived from this analysis (21°± 20°) is quite uncertain due to the 
nearly face-on orientation of NGC 4314, however it is consistent with the value found in PI 
(23°± 8°), assuming that the outer isophotes obtained from digitally stacking the POSS O 
and E plates are intrinsically circular. However, the outer isophotes may be intrinsically 
oval, since outer rings and pseudorings located at the OLR are generally elongated (Buta 
1993). An argument based on the TuUy-Fisher relationship (see Kenney et al. 1993) gives 
27° ± 4°. We choose z = 21° to determine the CO-derived rotation curve. 

In an effort to extend the rotation curve outward from the OVRO limit, we have 
adopted the rotation curve of Combes et al. (1992) for r > 10" . This curve is obtained from 
a mass model of the galaxy based on a K-band image, which presumably accurately traces 
the stellar mass in the galaxy. We choose the Combes model derived rotation curve over 
that from Quillen et al. (1994), since the agreement with the CO velocities in the center of 
the galaxy is better. 
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Fig. 13a shows our OVRO and the Combes et al. (1992) rotation curves. Agreement 
between the two rotation curves is good for r < 7" . The CO curve flattens at r = 7" with 
Vmax ~ 175 km s~^. Assuming that the ring is co-planar with the plane of the galaxy, has a 
radius of 6'/7, and a circular velocity of 175 km s~^, the period of rotation for the gas in the 
ring is 11 My. 



4.2. Deviations from Circular Velocity 

Deviations from circular velocity can be seen directly on the velocity map shown in 
Fig. 12. There are regions associated with the dust lanes just outside the nuclear ring (e.g., 
Fig. 5) where the isovelocity contours bend sharply by ~ 90°. To study these regions in 
more detail, wc obtain a velocity residual map by producing a model galaxy velocity field 
from the rotation curve. Fig. 13a, and subtracting it from the observed CO velocity field. 
This is shown in Fig. 14, overlaid on the CO intensity map. 

This comparison shows that the largest deviations from circular motion occur in two 
symmetrically located spurs, just outside the ring and near where the dust lanes merge with 
the ring. These spurs are located along the minor axis, where, within a disk, non-circular 
motions imply radial motions. If the dust lanes are truly along the leading edge of the bar, 
and the spiral arms are trailing, then the southwest side of the galaxy is closest to us and 
these radial deviations imply inflow. These local inflow speeds reach their maximum of 90 
km at r ~ 10" and gradually decrease as the nuclear ring is approached and crossed. 
Inflow along the leading edge of a bar is predicted by theory and models (e.g., Roberts et 
al. 1979, Athanassoula 1992b) but is rarely seen as clearly as it is in NGC 4314. 

We can estimate the mass inflow rate into the ring. The CO flux in the two inflow 
spurs (±7" from the nucleus) is 50 Jy km s"^, which corresponds to a gas mass of M(H2 -|- 
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He) = 7 X 10^ Mq. Adopting a mean inflow speed of 30 km s ^ implies that the gas now at 
7" will reach the nucleus in 336 pc/30 km s"-*^ = 10'^ yr. Thus the average inflow rate over 
this time is 7 x 10^ Mq/W yr = 7 Mq yr'^ 

Strictly speaking, the inflow speeds measured in Figure 14 are local streaming motions, 
and do not necessarily represent net inflow speeds. Gas along bars has the largest inward 
motion where it piles up and is easiest to detect, but is predicted to move outwards on 
other parts of its orbit, where the gas surface density is lower (Athanassoula 1992a). This 
is a greater concern far out along bars than just outside the ILR ring, and we believe that 
the net inflow speeds are close to the local streaming motions in this part of NGC 4314. 

From the thermal radio continuum flux, G-B estimate an extinction-corrected Ha flux 
for NGC 4314 of 1.5 x 10~^^ erg s~^ cm~^, which corresponds to a star formation rate of 
0.16 Mq yr~^, using an extended Miller- Scalo IMF as in Kennicutt (1983). Since this rate 
is much lower than the inflow rate, the gas mass in the ring may be increasing with time. 
And, since little CO is seen in the dust lanes, this process cannot continue much longer. 

In Fig. 15b the velocity deviation map is compared to the Pogge (1989) Ha image. 
Luminous Hll complexes are located within 1" of the terminus of each inflow spur. 
This suggests that gas flowing inward along the bar collides with gas already in the 
ring, triggering star formation. The Hll regions appear upstream from the inflow zones, 
suggesting that it is gas already in the ring that is involved with star formation. These 
regions may be similar to the bar end in M83, where two gas streams meet and one of them 
undergoes star formation (Kenney & Lord 1991). 

The existence of CO near the nucleus implies the need for a collection mechanism. In 
P2 we presented evidence for the existence of an 8" long nuclear bar, oriented parallel to the 
primary stellar bar. If this nuclear bar is a real bar-like dynamical feature, then gas flow 
patterns of the type thought to be associated with bars should be present. Fig. 16 shows 
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that deviations from circular motion (the same differential velocities presented in Fig. 14, 
but with AV = 5 km s"-*^) are seen along the nuclear bar. However, the velocity pattern is 
not everywhere consistent with the streaming motions generally associated with bars. The 
signal-to-noise ratio in this central region, where CO emission is weak, may be too low to 
show the streaming motions clearly. 



4.3. Resonances 



From the rotation curve in Fig. 13a, we produce standard resonance curves (c.f. 
Devereux et al. 1992), and plot them in Fig. 13b. Because the angular velocity of this 
bar is unknown, we assume a pattern speed of flp = 3.5 km s~^ arcsec"^, which places 
corotation at the end of the primary bar. A /ie-i map of NGC 4314 (Fig. 16) shows that 
it has an outer elliptical distribution of newer stars with r = 13" oriented perpendicular 
to the primary stellar bar (PI). NGC 4314 also has a nuclear bar with a 4" radius (P2) 
oriented parallel to the primary stellar bar. These orthogonal structures suggest nested 
resonances (Combes 1988). The resonance curves in Fig. 13b would predict OILR at r = 13 
and IILR at r = 4" for a bar pattern speed 0,p — 3.5 km s~^ arcsec"^, corresponding to 72 
km s~^ kpc~^. This places the outer Lindblad resonance (OLR) at r ~ 98", coincident with 
the elliptical outer isophotes discussed in PI. These elliptical isophotes have their major 
axis oriented perpendicular to the primary stellar bar, as do some rings and pseudo-rings in 
other galaxies that are located at the OLR (Buta & Crocker 1993). 



5. Comparisons with other Galztxies 
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A concise review placing NGC 4314 in the context of other barred galaxies with similar 
observations can be found in Kenney (1995). Here we offer a few specific comparisons. 
In the central regions of MlOl (Kenney et al. 1992), NGC 6951 (Kenney et al. 1992), 
NGC 1530 (Regan, Vogel, & Teuben 1995), and MlOO (Rand 1995; Sakamoto et al. 1995), 
CO maps show spiral-shaped features extending out from a circumnuclear ring or partial 
ring of Hll regions. In many of these examples the CO is strongest at the intersection 
of the ring with the dust lanes along the leading edge of the bar. This morphology is 
particularly striking in NGC 3351, where the CO is concentrated into "twin peaks" located 
symmetrically about the nucleus oriented perpendicular to the large-scale bar (Kenney et 
al. 1992). These peaks have been interpreted as regions of orbit crowding, where inflowing 
gas from the bar meets gas on more circular orbits near the ILR (Kenney et al. 1992). 

Our data for NGC 4314 show spiral-shaped inflow spurs, but they are weaker than 
in most of these other galaxies. NGC 4314 also does not have a "twin peak" morphology. 
There are modest CO peaks located where the inflow spurs meet the ring, but the strongest 
CO peaks are located elsewhere in the ring. The variety in gas morphology between galaxies 
is probably due in part to dynamical differences, including variations in the degree of central 
concentration and the strength of the bar. The rotation curve and the pattern speed of the 
stellar bar determine the number, location, and separation of the ILRs. Gas behavior also 
depends on whether there are small scale nuclear bars (Shlosman et al. 1989; Wozniak et al. 
1995), and the relative position angles of the bars, if they are independently rotating. Some 
of the variety may also be due to evolutionary differences. Timescales are relatively short 
in circumnuclear regions; starburst timescales are ~ 10^ — 10* yr, and dynamical timescales 
are ~ 10^ yr. The relative weakness of the CO in the inflow spurs of NGC 4314 compared 
to these other galaxies may mean that in NGC 4314 the gas in these dust lanes has already 
been driven into the center. The observation that the molecular gas in NGC 4314 lies 
inside the star forming ring lead Combes et al. (1992) to suggest that the ring is shrinking 
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with time due to dynamical friction from the bulge stars. As noted previously, there is an 
elliptical distribution of newer, but non-ionizing stars outside of the ring of H ll regions 
(PI). There is a color progression along this ellipse, such that stars closer to the ring are 
bluer (PI). Since our OVRO data show very httle CO along this eUipse, we conclude that 
this color variation is due to age rather than extinction. Therefore, younger stars are found 
closer to the Hll region ring, supporting the idea of a shrinking ring of star formation. 
This observation, along with the observation that essentially all of the molecular gas in the 
galaxy is found in or inside the ring and the fact that NGC 4314 is extremely deficient in 
atomic gas (G-B), lead Combes et al. (1992) to hypothesize that NGC 4314 is in a late 
stage of evolution, and that accretion to the NGC 4314 ring has stopped. Our observations 
show that inflow is in fact still occuring. However, since little CO is seen in the dust lanes, 
this process cannot continue much longer. 



6. Summciry 

The major observational results of this paper include 

1. We confirm the ring-like CO morphology first reported by Combes et al. (1992). 

2. We discover significant CO near the center of NGC 4314 without associated recent 
star formation. 

3. CO associated with dust along the leading edge of the primary stellar bar is detected 
at a 3-(T level. This association is confirmed for the dust bowl region by a strong 
correlation of CO intensity with extinction estimated from optical color maps. 

4. Much of the CO in and near the ring is coincident with dust. 
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5. We identify the nested, orthogonal resonances in NGC 4314. We identify the IILR 
with the end of the nuclear bar, OILR with the outer extent of an elliptical distribution 
of less recent star formation, corotation near the end of the large-scale bar, and 
the OLR with the eUiptical outer isophotes oriented perpendicular to the bar. The 
strongest star formation occurs between the IILR and OILR and closer to the IILR. 

6. A color gradient, from blue near the IILR to red near the OILR, is likely due to stellar 
age differences, not reddening. CO emission is strong near the IILR and weak near 
the OILR, hence the expected reddening gradient is the opposite of the observed color 
gradient. 

7. Radial inflow at speeds of 20-90 km s~^ is detected in two spurs of molecular gas 
located just outside the star-forming ring. Significant Ha emission occurs where the 
inflow spurs meet the ring. 

Benedict thanks Jeff Achtermann, Antonio Garcia-Barreto, Jim Higdon, and Hong Bae 
Ann for useful discussions. Melody Brayton provided invaluable assistance with LaTex. 
We thank Rick Pogge for providing an Ha map of NGC 4314. Benedict acknowledges 
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Figure Captions 



Fig. 1. — An unsharp masked image of NGC 4314 from PI. This high-pass map emphasizes 
regions within the galaxy where luminosity gradients are changing rapidly. It brings into 
prominence dust lanes (white) and (for r ~ 6") confirmed regions of recent star formation 
(black). The sharp gradient of the nucleus also shows up as black. Note the 26" diameter 
oval-shaped distribution at p. a. ~ 58°, and the dusty region (the "dust bowl") along the bar, 
southeast of the nuclear ring at RA = 12^22™32!7, Dec = +29° 53' 27" (2000). 

Fig. 2. — The uniformly weighted channel maps. The velocity interval is 13 km s~^. Channel 
9 corresponds to V — 1075.5 km s~^, while channel 24 is V = 880.5 km s~^. The contour 
levels are 3, 4, ... ,10 times our uniform map noise level of 28 mJy beam~^. The dynamical 
center is marked with a cross in each channel. 

Fig. 3. — The uniform weight CO contours superposed on the gray scale CO map. The beam 
size is 2"3 x 2'/2. The contour levels are 0.1, 0.2, ... , 1.0 times the peak intensity of 1.3 Jy 
beam~^ km s~^. The Sa level is 0.153 Jy beam"-*^ km s~^, shghtly higher than the lowest 
plotted contour. 
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Fig. 4. — The natural weight CO contours superposed on the gray scale CO map. The beam 
size is 3'.'1 x 2'.' 7. The contour levels are 0.1, 0.2, ... , 1.0 times the peak intensity of 2.74 
Jy beam~^ km s~^. The 3a level is 0.165 Jy beam"^ km s~^, slightly lower than the lowest 
plotted contour. The extended distribution lies along the primary stellar bar at p. a. = 148°. 

Fig. 5. — a) A uniform weight CO intensity contour map superposed on the central region of 
Fig. 1. The second lowest contour level is 5a above the noise level. The gray scale is encoded 
such that local excess (newer stars) is black and local deficit (dust) is white. Note that the 
elliptical locus of newer stars oriented along p. a. = 58° is devoid of CO. b) A natural weight 
CO contour map (Fig. 4) superposed on Fig. 1. Note coincidence of CO with the leading 
edges of the primary stellar bar and with the dust bowl. 

Fig. 6. — The uniform weight CO contour map from Fig. 3 superposed on an Ha + [N II] gray 
scale map (Pogge 1989). The gray scale is encoded to represent the strongest Ha emission 
as black. Map registration was determined by identifying the Wakamatsu & Nishida (1980) 
Ha-bright knots A and B. 

Fig. 7. — Detailed profiles along p. a. = 148°, the major axis of the primary bar. Surface 
magnitude and color indices {^b, fJ^B-n, and are from PI. The surface density of CO is 
from the uniform weighted map (Fig. 3). The transformation of CO to molecular gas surface 
densities is described in the text. We model the gas surface density peak to the SE with a 
gaussian. The overlay shows the range and quality of the fit. //Ha is in arbitrary units. In 
the bottom panel, the ratio Ha/dHj provides a qualitative and relative indication of the local 
SFR. The missing data in and //b-h near the nucleus are due to point spread function 
differences between the visible and short wavelength infrared (SWIR) bandpasses (see PI). 

Fig. 8. — Same as Fig. 7, except along p. a. = 58°, the minor axis of the primary bar. 

Fig. 9. — Same as Fig. 7, except along p. a. = 0°, through the HII region, knot A. 
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Fig. 10. — Profiles through the dust bowl perpendicular to the primary bar. The /xb, A*b-v 
and //v-i data are from PI. There are no J, H or K data for this region. 

Fig. 11. — Visual absorption (Ay) and CO surface density {Mq pc^^) are correlated for the 
gas and dust at the dust bowl. 

Fig. 12. — The naturally weighted velocity field superposed on a uniformly weighted gray 
scale map. The contour interval is AV — 20 km s~^. We have labeled several contours. 
The broken circle contour to the NW is = 1040 km s~^. The similar contour to the SE is 
V = 940 km s-\ 

Fig. 13.— a) The rotation curve from OVRO and Combes et al. (1992) data. The OVRO 
data extend out only to 9". The Combes data extend the rotation curve to r > 120". The 
errors in V arc 1 — a. The error bars in radius represent the beam size, b) The resonance 
curves. A pattern speed of 3.5 km s~^ arcsec"^ places the IILR at the end of the nuclear bar 
(P2) and the OILR coincident with the outer eUipse of newer star formation discussed in PI 
and shown in Fig. 14. 

Fig. 14. — Velocity residuals produced by subtracting the Fig. 13 rotation curve from the 
observed natural weight velocity field in Fig. 12. Residuals are superposed upon the Fig. 3 
uniform weight CO intensity gray scale map. The velocity interval is 10 km s~^. 

Fig. 15. — a) The natural weight CO velocity residual map superposed on the optical unsharp 
masked data, b) The natural weight CO velocity residual map superposed on Hq;+[N II] from 
Pogge (1989). The gray scale is encoded to represent the strongest Hq;+[N II] emission as 
black. Note that the deepest incursions of the decelerating gas terminate in close proximity 
to large H ll regions. 
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Fig. 16. — Surface I-band intensity (greyscale) from HST (P2) and CO velocity deviations 
from purely circular motion (contours). AV — 5 km s"^ . The nuclear bar is aligned along 
p.a= 148°. 



Fig. 17. — B-I surface color (//r-i) map of NGC 4314. Encoded such that darker to hghter 
represents bluer to redder. Note the discrete steps in the color progression along the feature 
identified with the OILR. 
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TABLE 1. Global Properties of NGC 4314 

TABLE 2. OVRO Interferometer Results for NGC 4314 



